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Generic spherically symmetric self-similar collapse results in strong naked-singularity formation. 
In this paper we are concerned with particle creation during a naked-singularity formation in spher- 
ically symmetric self-similar collapse without specifying the collapsing matter. In the generic case, 
the power of particle emission is found to be proportional to the inverse square of the remaining 
time to the Cauchy horizon (CH). The constant of proportion can be arbitrarily large in the limit 
to marginally naked singularity. Therefore, the unbounded power is especially striking in the case 
that an event horizon is very close to the CH because the emitted energy can be arbitrarily large 
. . . in spite of a cutoff expected from quantum gravity. Above results suggest the instability of the 

I CH in spherically symmetric self-similar spacetime from quantum field theory and seem to support 

. the existence of a semiclassical cosmic censor. The divergence of redshifts and blueshifts of emitted 

' particles is found to cause the divergence of power to positive or negative infinity, depending on the 

coupling manner of scalar fields to gravity. On the other hand, it is found that there is a special 
^ , class of self-similar spacetimes in which the semiclassical instability of the CH is not efficient. The 

■ analyses in this paper are based on the geometric optics approximation, which is justified in two 

' dimensions but needs justification in four dimensions. 

^ . PACS numbers: 04.20.Dw, 04.62. -fv 
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' The cosmic censorship hypothesis asserts that gravitational collapse of physically reasonable matter with generic 
initial condition does not result in the formation of naked sing ularity (NS) Q. More precisely, there are two versions 
of this hypothesis. The weak hypothesis states that all singularities in gravitational collapse are hidden within black 
holes. This implies the future predictability of the spacetime outside the event horizon. The strong hypothesis asserts 
that no singularity visible to any observer can exist. It states that all physically reasonable spacetimes are globally 
hyperbolic. Although the hypothesis has been an object of study over the last few decades, there is little agreement as 
to the final end of a general relativistic gravitational collapse. There is not enough evidence to prove the hypothesis. 
On the contrary, some solutions of the Einstein equation with regular initial conditions evolving into spacetimes which 
contain naked singularities have been found. 

Up to now, it has been reported that such semiclassical effects as particle creation by naked-singularity formation, 
which are analogous to the Hawking radiation 0, might prevent the formation of NS. The first researchers to give 
much attention to this possibility were Ford and Parker. They considered a globally naked-singular spacetime, in 
' which the weak hypothesis is broken, and calculated quantum emission from a shell-crossing singularity to obtain a 
finite amount of flux 0]. On the other hand, Hiscock, Williams, and Eardley obtained diverging flux from a shell- 
focusing NS which results from a self-similar implosion of null dust Q . Subsequently, these semiclassical phenomena 
have been investigated in the collapse of the self-similar dust 0, Q , non-self-similar analytic dust , and self-similar 
null-dust fluids Q . Typically, the power of quantum emission from shell-focusing NSs diverges as the Cauchy horizon 
(CH) is approached. For a recent review of the quantum/classical emission during naked-singularity formation, see j^. 
In addition, such an explosive radiation by naked-singularity formation can be a candidate for a source of the ultra 
high energy cosmic rays or a central engine of 7-ray burst [ij . 

As the examples given above show, it is known that generic spherically symmetric self-similar collapse results in 
strong naked-singularity formation jlll Il2l |. Among such self-similar models, the general relativistic Larson-Penston 
(GRLP) solution would be one of the most serious counterexamples against the cosmic censorship hypothesis in 
the sense that the existence of pressure is taken into account 0,^3- Moreover, the convergence of more general 



spherically symmetric collapse to the GRLP solution have been reported both numerically and analytically [15j as a 



een rep c 

realization of the self-similarity hypothesis proposed by Carr [1^, [13 • The discovery of the black hole critical behavior 
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also shed light on a self-similar solution as a critical solution (see jig for a review). We can say with fair certainty 
that self-similar solutions play important roles near spacetime singularities. Several studies have been done resulting 
in a complete classification of self-similar solutions so far (see for a review). 

Motivated by the above, this paper is intended as the investigation of particle creation during the naked-singularity 
formation in self-similar collapse and of the resulting instability of a CH. It is shown that irrespective of the details 
of the model, a diverging energy flux is emitted from a naked shell-focusing singularity forming in generic spherically 
symmetric self-similar spacetime. The power and energy of particle creation are calculated on the assumption that the 
curvature around the singularity causes particle creation, which was proven by |l9j | , at least in the case of self-similar 
and analytic non-self-similar dust models. Because collapsing matter is not specified, the results can be applied to 
several known models of self-similar collapse. Our analysis is regarded as a semiclassical counterpart of l23|, in which 
the stability of the CH in self-similar collapse was tested by a classical field. 

This paper is organized as follows. In Sec. a class of self-similar spacetimes is introduced and some features of 
the class of spacetimes to be used in later discussion are extracted. In Sec. IIIII null geodesic equation is considered to 
calculate the "local map" of null rays, which plays a central role in the estimation of the power of emissions. In SeclTvl 
the power and energy of the created particles are estimated, while redshifts are the focus of Sec. ^ Then several 
underlying relations among the local map, redshift, and power are found in Sec. I VII General results obtained in the 
preceding sections will be applied to the various models of collapse in Sec. IVIll Sec. IVIlJ is devoted to a summary 
and a discussion. In an Appendix, the same results arc derived in terms of another coordinate system, which is useful 
in applications. The signature (— , -I-) and units in which c — G — h — 1 are used below. 



II. SPHERICALLY SYMMETRIC SELF-SIMILAR SPACE-TIMES ADMITTING A NAKED 

SINGULARITY 



In this article a class of spacetimes which are spherically symmetric and admitting a homothctic Killing vector 
field ^, which satisfies £^gfj,i, = 'i.g^v, is considered. The line element of this class of spacetime in an advanced null 
coordinate system is written as 



(2.1) 



where x = v/R, dfl^ is the line element of a unit two dimensional sphere, and the homothetic Killing vector field is 
of the form ^ = vdy + Rdn. In this spacetime, the geodesic equation for an outgoing null ray is written as 



dv 
dR 



'^gvi 



xf{x), 



where 



Equation H2.2|l can be written also as 



xgvv 



dx 
dR 
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which is integrated to give 
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(2.2) 



(2.3) 



(2.4) 



(2.5) 



where xq and Rq are constants which are related as Rq — R(x — a;o). The constant xq is chosen as xq < x^ and 
xoj^O. 

What we have to do first is to extract features of f{x), which determines the spacetime structure. The Misner-Sharp 
mass in this spacetime is given by 



m(«,i?) = :^(l-V^i?V^i?) = ^ 



1 



x^Pgv 



The regularity of the center R = in the region v < and the absence of a trapped or a marginally trapped surface 
for < R and v < are assumed. The latter condition is W^RV^R > for all x e (— oo,0], which is written as 
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g^y < for all x € (— oo,0] in the present case. The inevitability of a curvature singularity at the origin v = R = 
can be shown except for a flat spacetime . In this article we consider self-similar spacetimes with a globally naked 
singularity, of which existence breaks the weak version of the cosmic censorship hypothesis. One of the possible causal 
structures of the naked-singular spacetimes is depicted in Fig. ^ The coordinate v is set to be the proper time along 
the regular center to remove a gauge freedom v V{v), so that lima;^_oo ffuw = —1- When m/R^ is required to be 
finite at the regular center, the function f{x) behaves as 

/ ~ 2/x a,s X ^ — OO. (2.6) 

The quantity m/R^ must be finite also in the limit v ^ {) for fixed R (> 0) so that 

/ 0{\xf) as X ^ 0, (2.7) 

where /3 < — 1. When there are positive roots of the algebraic equation f{x) = 1, it can be shown that the curve 
X = a;+ is a CH, as we will see in Appendix ^ where x'^ is the smallest root. The differentiability of the metric 
fimction / is assumed to be as follows: 

J e C" ((-c5o,a;+)) , /eC^- &i x ^ x+ . (2.8) 

The former condition guarantees the existence and uniqueness of geodesies in this system. It is also assumed 0| that 

f'{x+) < 0. (2.9) 
The schematic plot of the function f{x) is shown in Fig. I^Ja). 



III. LOCAL MAP 



To estimate the power of particle creation just before the singularity occurs, the pole at a; = should be extracted 
from the integrand in Eq. H2.5|l as follows: 
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where 



7 = x+f'{x+), 
F*{x) = F{x) 



(3.1) 



(3.2) 



"f{x — a;+) 



The constant i?o in Eq. I|2.5|) . which parameterizes solutions, is related to the time Vc = v{R — Q) < when the 
outgoing null ray emanates from the regular center as follows: 
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Taking the limit of i? ^ (w < 0), following relation is obtained: 



\xo\r 



I = exp 



f{x' 



X' [fix') - 1) 



dx' 



(3.3) 



where condition (|2.6|l ensures the convergence of the integral in Eq. (|3.3|l . 
Combination of Eqs. H3.1|) and H3.3|l yields 



R = C{R, x) {v+ {R) - v) ^''^ Ue, 



(3.4) 
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where 



v+{R) = x+R, 
C{R,x) = -IxoT^/"^ [(a;+ - a;o)-R]"^^'^exp 



F*{x')dx 



Before turning to the derivation of the local map, a few remarks should be made concerning the function C{R,x). 
Due to condition H2.8|l . C{R,x) converges to some finite constant in the limit cc — > a;+ for fixed R. The dependence 
of C (i?, a;) on a;o is only an apparent one as 



dC 
dxQ 



0, 



(3.5) 



which we use in Sec. I VII 

Now let us consider a pair of ingoing and outgoing null rays such that the latter is the reflection of the former at 
the regular center. An observer who rests at i? = 9^ will encounter the null ray twice so that we denote the time of 
first encounter by vi and that of the second by V2. The relation between vi and W2, which we call the local map, is 
obtained from Eq. 1)3.4(1 : 



Vl = 



where 



c{m,x2: 



X2 



{v+{^)-v2y 



V2m, 
-1/7. 



(3.6) 



(3.7) 



The time intervals —vi (> 0) and t;+($H) — V2 are depicted schematically in Fig.[5fb). It is noted that the result does 
not change even if we choose a small value of 9\. Namely, the nature of the local map is determined by the behavior 
of null rays near the singularity. 



IV. POWER AND ENERGY 



A global map V — G{U) is defined to be a relation between the moments when one null ray leaves I~ and 
terminates at X+ after passing through the regular center (see Fig.P). Assuming the geometric optics approximation, 
one can obtain the power of emission as the vacuum expectation value of a stress-energy tensor by the point-splitting 
regularization from the global map , 



P = 



for a minimally coupled scalar field, and 



247r 



P 



3 fG^'^ ' 
2\G' 



1 rc^ 

487r I G' 



'G' 



for a conformally coupled scalar field. 

Spherically symmetric self-similar spacetimes are not asymptotically flat in general. They therefore should be 
matched with an outer asymptotically flat region via a proper non-self-similar region. This matching procedure is 
quite straightforward for dust collapse Although it seems to be necessary to solve null geodesic equations in 
such a "patched-up" spacetime, the main properties of the global map must be determined by the behavior of null 
rays passing near the point where the singularity occurs. This expectation has been confirmed in |T9l |. at least for 
the self-similar and analytic dust models. Therefore, one can safely assume that the global map inherits the main 
properties of the local map, such as the value of the exponent and differentiability. This means that from Eq. 
the asymptotic form of the global map would take the form 



G{U) = Vo-iUo-UfG,{U), 
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where the null rays U = Uq and V = Vq are the CH and the ingoing null ray that terminates at the NS, respectively. 
The function G^(U) is a regular function which does not vanish at the CH and a is the exponent of the local map, 
ai in Eq. 

In the case of a = 1, the leading contribution to the power of particle creation is calculated as 

_ 2G:\Uo) - G.{Uo)G':{Uo) 

SirGliUo) ' ^ ' 



12n \G4U0) 

so that the power remains finite at the CH. Unfortunately G^{Uq) and its derivatives cannot be known until the null 
geodesic equation is solved globally, so that one could not know the power of emission from only the information 
contained in the local map. In terms of redshift, a = 1 corresponds to the case that the redshift of a particle remains 
finite at the CH, as we will see in Sec. |V] 

On the other hand, in the case of a 7^ 1 the leading contribution is obtained as 

a' - 1 

P= 2' 4.3 

48^ {Uq -Uf 

for a minimally coupled scalar field. For a conformally coupled one, the power of emission is obtained by replacing the 
factor (a^ — 1) in Eq. (|4.3|) for (q — 1)^. The power is proportional to the inverse square of the remaining time to the 
CH. If a > 1, the power diverges to positive infinity for both minimally and conformally coupled scalar fields, while 
if < a < 1, the power diverges to negative and positive infinity for minimally and conformally coupled scalar fields, 
respectively. In terms of the redshift of particles, the case that a > 1 (0 < a < 1) corresponds to infinite redshift 
(blueshift) at the CH, as we will see in Sec. El The emitted energy can be estimated as 

a' -I 

P(U')dU' ^ -. (4.4) 
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Although the emitted energy diverges when the CH is approached, this divergence needs to be regarded carefully. 
The semiclassical approximation would cease to be valid when the curvature radius at some spacetime point inside 
star reaches the Planck scale. Here we make a natural assumption that such a situation happens at the center of a 
star dX V = — iqc j22| |. In the case of a > 1, it can be expected that for a ray emanating from the center at t; = — igc, 
the time difference Uq — U would be greater than the order of Iqq due to redshift, i.e., AU = Uq — U > tqc- Then 
energy emitted by the time Uq — AU is 

— 1 — 1 , ^ 

where -Bqg = I/^qg- If the factor (a^ — l) /(487r) is on the order of unity, the total radiated energy within the 
semiclassical phase is less than the order of Eqq, which would be of course much less than the mass of ordinary 
astrophysical stars. It would be better to say that a collapsing star would enter the phase of quantum gravity with 
most of its mass intact. Therefore, one could not predict whether a star which collapses to a NS evaporates away 
or ceases to radiate at its final epoch. This has been pointed out in after careful investigation. We should not 
overlook that this feature is much different from that of black hole evaporation, in which quantum gravitational effects 
appear after a black hole loses almost all its mass. 

One can recognize, however, that if a 3> 1 the radiated energ y co uld be large. This situation is realized in the limit 
to marginally NS, in which the CH and event horizon coincide |24| . To illustrate the unbounded increase of a in this 
limit, we have to look deeper into the causal structure, which is determined by the function /. We order the positive 
roots of the equation f{x) = 1 as < xi = 2:+ < .X2 < ■ ■ ■ < Xn, where we count multiple roots as one root. The 
existence of Xa satisfying Imix^xa fi^) = with a;„ < Xa and the continuity of / in the region xi < x < Xa are 
assumed. In the region x S (a;„, Xa), dR/dv > along the null geodesies and Ymix^xc, dR/dv — from Eq. H2.2|l . This 
implies that outgoing null rays in this region are to turn back in the direction of the singularity eA, x = Xa and that 
the curve x = a;„ is the last outgoing null ray which can escape to infinity. That is to say, x = Xa and 
are the apparent horizon and event horizon, respectively. Hereafter, we shall concentrate on the case of n = 2. The 
function f{x) would be written as 

f{x)-l = Mx){x-X+){Xe-Xr , 0<X<Xa, (4.6) 
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where f*{x) is a function which satisfies f*{x~^) < and m is some positive integer. The exponent of the factor 
(x — x+) is restricted to unity because of the condition f'{x~^) < and the differentiabihty / € at a; = a;+. With 
Eq. 1)4. 6|l . the exponent of the local map is calculated as 

^ ^x+f'ix+) = 'x+nix+) ix, - x+r ' ^^■^^ 

to show that ai can be arbitrarily large in the limit x^ —^Xe- 



V. REDSHIFT 



The estimation of redshift of the radial null ray would help us understand the behavior of the power and would 
be necessary for discussing the validity of geometric optics and semiclassical approximations. Hereafter the tangent 
vector of the null ray is denoted by fc^ = dx^/dA, where A is an afhne parameter. 

With the null condition fc^fc^ = 0, the w-component of the geodesic equation fc^'V^fc'^ ~ leads to 



dfc" (k'^y f 1 dgyR X 1 dg. 



dX 



+ t:- 



R \gvR dx 2 gyfi dx 



0. 



Furthermore by using the relation 

k'"{x) is integrated to give 

fc"(a;) 



dX^ ^ 



h'" 



exp 



where 



F{x')dx' 



Fix) 



exp 



1 dg^j 



xgvv 
'^9vrJ dx 



F*{x')dx' 
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-(l+7)/7 



F*{x) EE F{x) + 



g„R dx ' 1 - f \x f dx 
1 + 7 1 



7 X ~ x^ 



The constant xq, which is set as xq < x^ , and the constant fcp are related as fcg — k'"{xo). 
The constant fcp is related to fc" = k'"{R~ 0) as 



exp 



F{x')dx' 



Combination of Eqs. (|5.1|l and (|5.2() yields 



where 



C{R,x) = [ix+ -io)R] 



(l+7)/7 



exp 



F*{x')dx' 



(5.1) 



(5.2) 



(5.3) 



Now, let us consider time-like observers who rest at i? = and R = U\ {dO = dcj) = 0). The observed frequency 
is given by cD = —u^k'^ — \/\gvv{x)\k'"{x), where is the four-velocity of observer. When cDi = \mix^-oo i^ix) and 
LU2 = \im^^^+ uj{x) are defined, Eq. (|5.3|l yields 



UJ2 



gvv {x^ 



.g™(-oo) 



Cim,X2)iv+iUl)-V2) 



ai — 1 



(5.4) 



Thus we see that if ai > 1 (0 < ai < 1) the redshift (blueshift) of emitted particle diverges at the CH, while it 
remains finite if ai = 1. The relation between the redshift derived above and the local map will be presented in the 
next section. 
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VI. RELATIONS AMONG THE LOCAL MAP, POWER, AND REDSHIFT 

There would be a relation between the local map and redshift because the local map describes a kind of time 
delay. Since the asymptotic behavior of the local map and redshift in the limit x ^ is considered here, the time 
dependence is omitted as C {R, x) ^ C (R) . From Eq. H3.5|) , the relation 

dv2 _ 7C(9^),^,+ ,^,_ (1+^)/^ 



dvi 9\ 



{v+im)-v2)^'+''>^-' (6.1) 



holds to give an alternative definition of redshift. Indeed, the time dependence in Eq. can be replaced with the 
ratio of fc" by Eq. H5.3|l as 



dV2 

dvi 



where we set xq = xq in the evaluation of the integral in C to derive Eq. I|6.2fl since C does not depend on xq from 
Eq. H3.5|l . Equation (|6.2|l can be written as 

= (6-3) 



where dr^ = ^J\gvv\dvi (i = 1, 2) is the proper time measured by the observer. This relates the time delay and redshift 
to reveal that the redshift essentially corresponds to the local map and also to confirm the consistency of the analyses 
in Sees. Uni and El 

There exists a plausible relation also between the power of emission and the redshift of particles as mentioned in 
Sec. II VI In the case of a > 1 (0 < a < 1), the power and redshift (blueshift) diverge at the CH from Eqs. H4.3|) 
and H5.4|l . while the power and redshift remain finite at the CH in the case of a = 1. We may, therefore, reasonably 
conclude that the divergence of the redshift or blueshift at the CH causes that of the power. 

VII. EXAMPLES 

In this section, we will take examples to illustrate how the features obtained in the previous sections are realized 
in concrete models. Since several models are written in the diagonal form of a metric tensor, the formulation and 
notations for the diagonal form of a metric tensor are developed in Appendix IbI 

A. Minkowski spacetime 

Although neither singularity nor horizon exists in Minkowski spacetime, one can test the formalism by applying it 
to this trivial spacetime. The line element is written as 

ds^ = -dv"^ + 2dvdR + R^dn^, 

or 

ds^ = ~dt^ + dr^ +r^dVP. 
From definitions (|2.3|) and (|B2|) . one obtains 

2 1 

f{x) = -, w±(z) = ±-. 
X z 

The roots of the algebraic equations f{x) = 1 and w±{z) — 1 are given as a;+ = 2 and z^ = ±1, respectively. It must 
be noted that the function f{x) satisfies conditions (|2.6I) - (|2.9|) . which have been assumed to derive the local map. 
According to Eqs. (|3.2I) . 13.7|l . (|B4I) . and HB8I) . one can easily check that the exponents of local maps, ai and a2, are 
unity. This fact just tells us that the redshift and power remain finite in a flat spacetime. 
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B. Vaidya Solution 

The Vaidya solution describes the coUapse of a null-dust fluid 25j. The global map for the self-similar Vaidya 
collapse was derived in 0- 

The line element in the self-similar Vaidya spacetime is written as 

ds^ = - (1 - m{x)) dv^ + 2dvdR + R^dft'^, 

where m(x) = for x < and m(x) = 2/ia; for a; > 0. The constant /x is restricted as < /i < 1/16 for the nakedness, 
so that the spacetime with /i = 1/16 corresponds to a marginally naked-singular one. The function ](x) is written as 

/(^) = 



x{\ — mix)) ' 

The roots of algebraic equation ](x) — 1 are given as 



1 - VI - 16Ai J 1 + VI - 16m 

and — 



4^ 4/x 
What has to be noticed is that / satisfies conditions (|2.6(l - (|2.9|l . The exponent is calculated as 



^ 1 _ 1 + VI- 

~ a;+/'(a;+) " 2V1 - 16^ ' 

This exponent coincides with that of the global map which was obtained in 8]. One can see that lim^^o ai = 1 and 
lim^^i/16 0L\ = oo, where the former and latter correspond to the limits of Minkowski and marginally naked-singular 
spacetime. This example is a good illustration of the efficiency of the local map and the divergence of a\ in the limit 
to marginally NS. 



C. Roberts Solution 

The Roberts solution describes the self-similar collapse of a niassless scalar field j2^. The line element is given as 

ds'^-ix- ^ ^(^)^'(^) \dv'+ , ^ dvdR + R'dn', (7.1) 
\^ ^l + h^x)J V^ + h^{x) 

where h(x) = for a; < and h{x) — ax for a; > 0, so that the region with negative v is flat. The constant a is 
restricted here as \a\ < 1/2 for the causal structure such that the spacetime has a time-like NS as in Fig. 13 The 
function f{x) is written as 



x(Vl + h^- 2hh') 



The conditions on f{x) are again satisfled. The algebraic equation f{x) = 1 has a positive root x+ = 2/Vl — 4(7^ so 
that the exponent of local map is calculated as 

for \a\ < 1/2. Thus we see that the Roberts solution provides a non-trivial example of spacetime in which the power 
of particle creation remains flnite at the CH in self-similar collapse. 



D. Lemaitre-Tolman-Bondi Solution 



The Lemaitre-Tolman-Bondi (LTB) solution describes the collapse of a dust fluid 27]. Although both global and 
local maps were derived in 5] and 13] respectively for self-similar LTB collapse, the exponent in 5] is reproduced 
from the formalism in Sec. IIIII 
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The line element of self-similar LTB spacetime (for example, see 0) is 



1 — az/3 
(l-az)i/3 



where the constant a is related to a "mass parameter" A as a = fv'A. The constant A is restricted to the range 
< A < 6(26 — 15-\/3) = Am, where the latter inequality is imposed by the nakedness of the singularity [S^, so that 
the spacetime with A — A„i corresponds to marginally naked-singular one. One obtains the function w±{z) according 
to definition (IB2II as 



w±{z) = ± 



1 — az/3 
z {I — az) 



1/3- 



(7.2) 



1/3 

The required conditions on w± in calculating the local map are satisfied. When a new variable y = {1 — az) 
introduced, Eq. H7.2() can be written as 



wl{z) = 1 



[g+{y)+g^{y)X 



where g±{y) = ?>y'^^ay'^ — ?>y^2a, so that the roots of algebraic equations ^±(2:) — 1 correspond to those of 5=p(y) — 0. 
Using the chain rule d/dz = {dy /dz)d/dy, one obtains 



2(l-3a|)ff^(a^) ^ 

3a±g±{a^) 



1/3 

where a± = (1 — az^) and the prime denotes differentiation with respect to the argument of the function. The 
exponent of local map is calculated as 



a2 



_ z w'_{z ) o?_g\{a^) 
z+k;^(z+) o?^g'_{ctJ) 



to coincide with the one derived in |a,UJ|- This exponent becomes close to unity as A — > and increases monotonically 
with A to diverge to infinity as A ^ Am- 



E. General Relativistic Larson-Penston Solution 



As the last example the general relativistic Larson-Penston (GRLP) solution, which describes the self-similar 
collapse of a perfect fluid IT^. is considered. For the present, it may be useful to review the GRLP solution 
and its importance, although we have mentioned them in the Introduction. The equation of state must have the 
form of P = kp from the requirement of self-similarity, where is a constant. The GRLP solution represents a 
naked-singularity formation in the range < fc < 0.0105, where the upper bound is imposed by the nakedness of the 
singularity. This solution is interesting because it provides the first example in which the pressure does not prevent 
the formation of a NS. Moreover, the convergence to the GRLP solution of more general solutions near the central 
region of stars have been strongly suggested numerically and supported by a mode analysis as a realization of the 
self-similarity hypothesis I?]. From the above, the GRLP solution will be a strongest known counterexample against 
the cosmic censorship hypothesis. 

Because the GRLP solution is a numerical one, an explicit expression of the exponent of local map could not 
obtained analytically, although it is unlikely that a is equal to unity. To show rigorously that the power of emission 
is proportional to the inverse square of the remaining time to the CH and that the constant of proportion diverges in 
the limit to marginally NS, the dependence of exponent a on k should be clarified numerically |29j. 



VIII. SUMMARY AND DISCUSSION 



We have been concerned with a quantum mechanical particle creation during the naked-singularity formation 
in spherically symmetric self-similar collapse. The power, energy, and redshift of emitted particles are analytically 
calculated on the assumption that the curvature around the singularity causes particle creation and the metric function 
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/ is around the CH. As a result, in the generic case in which the exponent of the local map a ^ 1, the power 
has been found to diverge as P oc (i^h — where (tc/i ~ i) is the remaining time until a distant observer would 
receive a first light ray from the NS. It is worth pointing out that the weaker differentiability of the CH leads to 
a different result, i.e., the power of emission has a different time dependence, although we have not looked deeper 
into such a possibility in this paper. The square inverse proportion of the power to the remaining time is due to the 
scale invariance of self-similar spacetimes. The constant of proportion has been found to be arbitrarily large in the 
limit to marginally NS. Therefore, this explosive radiation is especially striking in the case that the event horizon is 
very close to the CH because the emitted energy can be arbitrarily large in spite of a cutoff expected from quantum 
gravity. We go on from this to the conclusion that if the back reaction to a gravitational field is taken into account, the 
semiclassical effect would cause the instability of the CH and might recover the cosmic censor in this limiting case. On 
the other hand, in the non-generic case in which a = 1 the power remains finite at the CH, so that the semiclassical 
instability of the CH seems not to be efficient for this special class of self-similar solutions. The collapse of a massless 
scalar field described by the Roberts solution indeed does correspond to this case. In addition, it has been found that 
the diverging redshift and blueshift cause the divergence of the power to positive or negative infinity, depending on 
the manner of the coupling of scalar fields to gravity. The divergence will be a criterion for the stability/instability 
of a CH in a gravitational collapse. Having suggested that the particle creation could cause the instability of the CH, 
we still have a long way to go before knowing the existence of a classical/semiclassical cosmic censor in more generic 
spacetimes. 
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APPENDIX A: FIRST OUTGOING NULL RAY FROM THE SINGULARITY 



We shall prove the absence of an outgoing radial null geodesic which emanates from the origin ti = i? = before 
the null geodesic x = . Let {v,R) be a point on a radial null geodesic I in the region < x = v/R < a;^ which 
emanates from v = R = 0. Then, one can say that 



dx 
dv 



1 



(1-J)>0. 



(Al) 



By the uniqueness of the solution of the null geodesic equation, I cannot cross x 



at points other than v = R = 0. 



Therefore above inequality says that x decreases as u — ^ +0 but is bounded from below by 0. Hence the limit 



lim x(v) 



exists and satisfies < x < x^. In the case of < a; < , 



R 



holds, but 



lim R{v] 



= hm — — 

J v^+O x(v) 



= 



(A2) 



(A3) 



lim — -— 



lim — -— 

v^+o R'{v) 



lim xf{x) 



xf{x) 



(A4) 



where the I'Hopital's rule is used in the second equality. This contradicts the assumption that a;+ is the smallest 
positive root of equation f{x) — 1. Next, we shall consider the case of a; = 0. The fact that the solution R converges 
to a finite nonzero constant as x — > from the condition (|2.7|l contradicts the assumption that I is a geodesic which 



emanates from v = R = Q. Thus the case of x 
null ray which emanates from the singularity. 



is also excluded. Thus we see that x = x^ is the first outgoing 
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APPENDIX B: LOCAL MAP AND REDSHIFT IN DIAGONAL COORDINATES 

The line element of the class of spacetimes discussed in Sec. I^in a diagonal coordinate system is written as 

(is2 = gu{z)dt^ + grr{z)dr^ + r^S^{z)dQ^, (Bl) 

where z = t/r and S" is a dimensionless metric function. The homothetic Killing vector field is of the form ^ — tdt+rdr- 
If one defines functions w±{z) as 



w±(z) 



1 

±- 

z 



9tt ' 



(B2) 



the roots of the algebraic equation 'w±{z) = 1 play important roles as do those of f{x) = 1 in Sec.^ The existence 
of a positive (negative) root of w-\-{z) — 1 (u)_(z) = 1) and the uniqueness of the root of w-{z) — 1 are assumed. 
When the negative root and the smallest positive root are denoted by z~ and 2+ respectively, the curves z = z^ and 
z = z^ can be shown to be the CH and the ingoing null ray that terminates at the NS. In addition, it is assumed 
that w'^i^z^) ^ 0. 

The null geodesic equations are integrated to give 



r 

— = exp 



W±{z')dz' 



exp 



z^ — z 



l/5± 



(B3) 



where 



W±{z) 
S± 

Wliz) = W±iz) 



z{w±{z) ~ 1)' 
S± = z^'w'j.{z^), 



(B4) 



1 



6±iz-z±)' 



and the signature + (— ) corresponds to outgoing (ingoing) null geodesic. The constants Zq and Tq are related as 



r{z = Zq) for an outgoing ray, while rg 
The constants are related to tc = t{r ~ 0) as 

±_ tc 

'''0 — I ±1 T ' 



r[z — Zq ) tor an ingoing one, where Zq is set as Zq < z^ and 7^ 0. 
w±(z')dz' 



J± — exp 



Combination of Eqs. (jB3|l and (jB5|l yields 



z'{w±{z') ~ 1) 



r = D±{r,z){t^{r)-tY"^U, 



(B5) 



(B6) 



where 



i±(r) 
D±{r,z) 



z^r, 



)r\ exp 



Wl{z')dz' 



Now, let us consider a pair of ingoing and outgoing null rays such that the latter is the reflection of the former at 
the regular center (r = 0,t < 0). An observer who rests at ?' = r will encounter the null ray twice, so that we denote 
the time of first encounter by ti and that of the second by t2. By Eq. ljB6|) . the ingoing and outgoing null rays are 
matched at the center to give the local map as 



D-iv, zi) 



5- 



where 



Zi 

a2 



U/v, (i = 1,2), 
6-/S+. 



(B7) 



(B8) 



12 



The redshift of a radial null ray is obtained in similar way. The t-component of equation fc^V^fc"^ = is integrated 
to give 



Lt,± 



exp 







W±{z')dz' 



exp 



Wl{z')dz' 



(B9) 



where 



W±{Z) = T TT- < (1 - 2W 1 ) ; 1 — 

2(l-w^^) V ^ ' gtt dz grr dz 

1 4- A , 1 

Wl{z) = W±{z) + 



-'± 

l + 5± 1_ 

5± z — z^ 



The constants z^ and fcp^ are related as k^^ — k^{zQ ) for an outgoing ray, while fcg = fc*(^o ) ^^''^ ingoing one. 



The constants /cq^ are related to A:* = k*{r = 0) as 



/cq^ — -T^, J± = exp 







W±{z')dz' 



Combination of Eqs. (|B9p and HB10|I yields 



k\z) = b^{r,z) (t^{r)-t)-^''-'^^"^ kl 



(BIO) 



(Bll) 



where 



1 r/ ± ~±N -ia+s±}/s± 



exp 



Wl{z')dz' 



Consider again the observer who rests at r = r and the pair of ingoing and outgoing null rays. The outgoing and 
ingoing null rays are matched at the center by Eqs. (|B7|I and l|Blip to give 



UJ2 



9tt{z ) 



D+iv, Z2) 



D+{r, Z2) 



D-ix, zi) 



{t+i.)-t2r--\ 



(B12) 



where lji = lim^^^^- \/\gttizi)\ fc*(zi) and 0)2 — lim^2_^2+ ^/\gtt(z2)\ k''^ {Z2) are the observed frequencies. 

There exists a plausible relation between the local map and redshift. From Eqs. ljB7(l and ljB12|l . one obtains 

dT2 _ LOl 

dri Qj2 ' 



where dr^ = \f\gu\dti (i — 1, 2) is the proper time of the observers. 
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FIG. 1: One of the possible causal structures considered in this article. A singularity occurs at the spacetime point (!7o, Vb) 
and is visible from X"*", where {U,V) are suitable double null coordinates. An outgoing null ray U = const, can be traced 
backward in time from to T~ , which turn out to be an ingoing null ray V = G{U). The outgoing null ray U — Uo and 
ingoing null ray V = Vq represent the CH and the null ray that terminates at the NS. 
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(a) (b) 

FIG. 2: (a) Schematic plots of f{x) defined by Eq. 12.31 for typical collapsing spacetimes which end in a NS or a black hole. 
Depending on the number of roots of f{x) — 1, which we denote by j, the causal structure of spacetime changes. The cases 
of j = 0, l,and 2 are depicted, (i) The case of j = 2 : f{x) (x > 0) is depicted by a solid line. The two roots are denoted 
by a;^ and Xe {x'^ < Xe)- The geodesies x — x^ and x — Xe represent the CH and event horizon, respectively. This kind of 
spacetime admits a NS. (ii) The case of j = 1 : f{x) {x > 0) is depicted by a dashed line. In this holds, i.e., the 

CH and event horizon coincide. This type of singularity is called marginally naked (MN). (iii) The case of j = : f{x) {x > 0) 
is depicted by a dot-dashed line. In this case the collapse ends in a black hole (BH). (b) A typical spacetime diagram of a 
collapsing body which ends in a naked singularity in {v,R) coordinates. A null ray which is reflected at the regular center 
and characteristic null rays in respective regions divided by horizons are depicted. The time intervals w"'"(5H) — «2 and ~vi in 
Eq. 1)8. 6|l are depicted. The dotted line is the world line of an observer at 7? = $R. 




FIG. 3: The conformal diagram of the Roberts solution Eq. 1)7.1(1 with \a\ < 1/2. The region « < is flat and the shading 
region is flUed with a collapsing massless scalar fleld. A time-like NS occurs at u = 0. 



